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ABSTRACT 

NGC604 is the largest Hn-region in M33, second only within the Local Group to 30 Dor, 
and is important as a laboratory for understanding how massive young stellar clusters interact 
with the surrounding interstellar medium. Here, we present deep (300ks) X-ray imagery of 
NGC604 obtained as part of the Chandra ACIS Survey of M33 (ChASeM33), which show highly 
structured X-ray emission covering ^70% of the full Ha extent of NGC604. The main bubbles 
and cavities in NGC604 arc filled with hot (kT — 0.5 keV) X-ray emitting gas and X-ray spectra 
extracted from these regions indicate that the gas is thermal. For the western part of NGC604 
we derive an X-ray gas mass of ^4300 Mq and an unabsorbed (0.35- 2.5 keV) X-ray luminosity 
of Lx = 9.3 x 10 35 ergs _1 . These values are consistent with a stellar mass loss bubble entirely 
powered by about 200 OB-stars. This result is remarkable because the standard bubble model 
tends to underpredict the luminosity of X-ray bright bubbles and usually requires additional 
heating from SNRs. Given a cluster age of ~ 3 Myr it is likely that the massive stars have not yet 
evolved into SNe. We detect two discrete spots of enhanced and harder X-ray emission, which we 
consider to be fingerprints from a reverse shock produced by a supersonic wind after it collided 
with the shell wall. In the eastern part of NGC604 the X-ray gas mass amounts to ~175OM . 
However, mass loss from young stars cannot account for the unabsorbed X-ray luminosity of 
Lx = 4.8 x 10 35 ergs _1 . Off-center SNRs could produce the additional luminosity. The bubbles 
in the east seem to be much older and were most likely formed and powered by young stars and 
SNe in the past. A similar dichotomy between east and west is seen in the optical, implying that 
a massive wall of neutral and ionized gas shields the dynamically quiescent east from the actively 
star forming west. 
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1. Introduction 

NGC604 is the second largest giant H n-region 
in the Local Group and has been studied through- 
out all wavelengths, from radio ( Churchwell & 



Goss Il999l: iTosaki et al. 



don et al. 



2007h. infrared CHig- 
2003 1. optical (iRosa k D'Odorico 



19821; iTenorio-Tagle et al.ll2000l). ultraviolet ( Rosa 



198ol: iKeel et al. Il2004lb to the X-rav regime 
Maiz-Apellaniz et al. 1 1200-41 ) . The main focus of 
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these studies was the investigation of the efFects 
that massive young stars impose on the surround- 
ing multi-phase interstellar medium. 

It is commonly believed that NGC604 is ion- 
ized by the radiation field and winds from a young 
and very massive stellar population. A detailed 
analysis of the stellar content indicates the Exis- 
tence of about 200 O and WR-stars with a WBL/O- 
ratio o f about 0.075 forissen et al. 1 19931; Hunter 
et al. Il996l: iGonzalez-Delgado k Perez I E OOOl ; 



Bruhweiler et al. 1120031 ). Although the st ellar 



content, the gas mass, and the average gas den- 
sity of NGC604 are different compared to other 
giant Hii- regions in, e.g., the LMC , IIZw4'), or 
NGC4214 (jTenorio-Tagle et al.ll20bll . the blister- 
like structures and overall Ha morphology result 
ing from stellar winds and/or SNe activity is 
similar among all these regions. In this re 
the similarity with 30 Dor in the LMC is the nost 
striki ng (see, e.g., Wang 19991 ; Townslev et 



2006). 



very 
gard 



al. 



Besides the warm ionized phase of the ISM in 
NGC604, the hot ionized medium (HIM) is an- 
other crucial, yet poorly investigated, phase of the 
ISM which tells us much about the interna en- 
ergetics, the wind-ISM interaction, the evolution 
of the bubble, and the sources which power it. A 
difficulty in the past originated from the fact that 
several X-ray bright superbubbles (logLx > 35) 
in the LMC show an excess X-ray luminosi y in 
comparison to predictions from the standard pub- 



ble m odel ( Castor et al. 1975; IWeaver et 



19771 ). There is general agreement that this ex- 



cess emission can be attributed to off-center SNRs 



interacting with the shell wa lls (|Chu fc Mac 
Il99d ; IWang fc Helfand |[l99lh . On the other bland, 

with 
onal 



X-ray dim bubbles are found to be consistent 
the standard model and do not require addit 



heating from SNRs ()Chu et al. 1119951 ). Another 



al. 



subdivision of stellar bubbles was introduce i by 
lOev I (1996) who distinguished between higm ve- 
locity bubbles, i.e., those whose expansion veloc- 
ity is underestimated by the standard modefl and 
those which are consistent with it. The emerging 
picture seems to be that X-ray bright, high vbIoc- 
ity bubbles are additionally heated by SNRs and 
that X-ray dim, slowly expanding bubbles do not 
require additional heating. 

With the current study, which is part o: the 
Chandra ACIS Survey of M33 (ChASeM33, Ipiu- 



2008), we present the first detailed 



cinsky et al 

X-ray analysis of NGC604 in order to explore 
the origin o: the X-ray emission and to test the 
wind-blown p ubble scenario. Altho ugh NGC604 
(L> = 817kpqV 



tect ed in previous maior X-rav surveys f Long et 



al Il996l; iMisanovic et al. 1120061; Plucinskv et 



Freedman et al. 1120011 ) has been de- 



al. 20081 ). th" ; only significant observation that has 



been publisl ed is a Chandra contour map of the 
central part ( Maiz-Apellaniz et al. 20041 ). which 
revealed tha; the bubbles seen in the optical are 
filled with X-ray emitting gas. However, the lim- 
ited sensitivity from this 90 ks observation (ObsID 
2023) inhibited further analysis. 

The papei is structured as follows: after an out- 
line of the data reduction steps for imaging (Sect. 

3troscopy (Sect. 2.2), we present the 



2.1) and spe 
main results 



derived from both techniques in Sect. 



3. In Sect. 4.1 we discuss the origin of the X-ray 
gas. By means of spectral modeling, we deter- 
mine electro i densities and filling factors of the 
HIM and derive X-ray gas masses for individual 
4.1.1). In Sects 4.1.2 and 4.1.3 
we cross-chei whether the stellar O/WR popu- 
lation in NGC604 can account for the observed 
X-ray gas m isses via continuous stellar mass loss 
served X-ray luminosity is consistent 
predicted from the " standard" model 



and if the of 
with the one 



of a stellar wind-blown bubble ( Weaver et al 



19771 ). Sect. 4 is concluded by evaluating whether 
SNe produce i a significant fraction of the observed 
X-ray gas nass (Sect. 4.2) and by presenting a 
ture of NGC604 which emerges from 
the multi-wE velength data analysis (Sect. 4.3). A 
summary is arovided in Sect. 5. 



2. Observations and Data Reduction 

This study is part of the Chandra ACIS Sur- 
vey of M33 (ChASeM33) and presents the first 
deep and high resolution X-ray images of the dif- 
fuse emissio i of NGC604. We chose 2 different 
ACIS-I fieldsL Field 2 from ChASeM33 (using Ob- 
sIDs 6378, 6B79, and 7402, -4' off-axis) and one 
archival field (ObsID 2023, with NGC604 at the 
aim point; s ee also dashed circ le in Fig. 1 of Plu- 
cinsky et al. (|2008J)). The total integration time 
was 300 ks. 

All ChASeM33 observations were performed 
with ACIS-j in 'Very Faint' (VFAINT) mode, 
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whereas the archival one was obtained with ACIS-I 
in 'FAINT' mode. All data were reprocessed with 
CIAO version 3.4.1.1 and CALDB version 3.3.0.1. 

The basic data reduction was carried out by fol- 
lowing standard procedures. First, a new level=2 
eventfile was created with acis_process_events 
by applying new background flags for the VFAINT 
and FAINT modes, correcting for charge trans- 
fer inefficiency (CTI) and time-dependent gain 
variations. Pixel randomization was removed 
to avoid degrading the resolution of the data. 
The resulting event lists were filtered for events 
with grade=0,2,3,4,6 and status=0, and the 
pipeline-generated good time intervals (GTI) 
were applied. No aspect solution correction was 
neededj. Positional uncertainties are < 1". 

Subsequently, we checked for background flares 
by extracting a background light curve from the 
events on the ACIS-I3 chip, after removing point 
sources detected by the task wavdetect. An itera- 
tive sigma-clipping algorithm was used to remove 
time intervals with count rates more than 3tr from 
the mean of each iteration, until all count rates 
were within ±3cr of the mean H The GTI-corrected 
total exposure time was 294 ks. 

2.1. Imaging 

We created X-ray images for each ObsID in 
the following diagnostic energy bands: 0.35- 
8.0keV (broad), 0.35-1. IkeV (soft), l.l-2.6keV 
(medium), and 2.6-8.0 keV (hard). Individual 
exposure maps for each ObsID and energy band 
were computed using the CIAO tools mkinstmap 
and mkexpmap. In order to ensure that the number 
of counts in each pixel of the exposure corrected 
image is proportional to the integrated number 
of counts of the incident source spectrum (over a 
given energy band), weighted instrument map^l 
were created assuming a spectrum representative 
of the diffuse X-ray emitting gas. We assumed an 
APEC thermal plas ma model with kT — 0.4 keV 
(|Long et al. Ill996h and an averaged column den- 
sity of TVtt = 1.1 x 10 21 cm~ 2 (Chur chwell & Goss 



19991 : iLebouteiller et al. 1120061 ). The script 



merge_all was used to produce a merged event 
list and merged exposure maps for the different 



energy bands. Merged images with integer counts 
(Fig. [T]) were created with dmcopy after applying 
bin and energy filters for the individual energy 
bands. Exposure-corrected images were created 
by dividing the (merged) images by their associ- 
ated (merged) exposure maps. All images were 
binned by 4 pixels to a resolution of ~2" to gain 
sensitivity and to compensate for varying PSFs at 
different off-axis angles. The X-ray images which 
were used for the contour maps (Fig. ^ and the 
RGB-i mage (Fig. [3t were a dditionally smoothed 
in DS9 (| Jove fc Mandel II2003T ) with a Gaussian ker- 
nel radius of 2 pixels, to help enhance the contrast 
between background and diffuse emission. 

For a multi-wavelength comparison of the dif- 
ferent phases of the ISM in NGC604, we also 
employed archival data collected at near infrared 
(NIR, Spitzer/TRAC, 8 fan), optical (HST/WFPC2, 
Ha, [Nil]), and far ultraviolet (FUV, GALEX, 
1500 A) wavelengths. Because we use the supple- 
mental images exclusively for morphological com- 
parison, standard data pipeline calibration was 
applied. The HST Ha image was also continuum 
subtracted using a scaled version of the V-band 
image. 

2.2. Spectroscopy 

Because we are mainly interested in the purely 
diffuse emission of NGC604, we removed point 
sources fr om every event li st using the ACIS 
Extract (jBroos et al.l 120021 ) r egions from the 



ChA SeM33 point source catalog (jPlucinskv et al 
2008). For proper background removal, we en- 



1 see http:/ /cxc. harvard. edu/ciao3.4/threads/arcsec_correction/ 

2 http: / / cxc. harvard, edu /ciao3 . 4/ threads / filter_ltcrv / 

3 http: // cxc. harvard, edu /ciao3 . 4/ threads / spectral, weights / 



sured that the background regions (boxes) were 
located at the same distance from the readout 
node as our source regions (to account for the row 
dependence of front-illuminated CCDs) . In order 
to improve photon statistics for the background 
data, we chose the sizes of these background re- 
gions to be approximately three times the size of 
the source region, yielding typically ~1200cts. 

The NGC604 pointings (Field 2 and the archival 
one) differ in pointing direction and roll. We there- 
fore did not produce a single merged event list and 
instead extracted source and background spectra 
from each individual event list. For this pur- 
pose we used specextract, which also creates cus- 
tom response matrix files (RMFs) and appropriate 
auxiliary response files (ARFs). The merging of 
RMFs and ARFs was done using the FT00LS tasks 
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addarf and addrmf and by applying a weighting 
factor based on the GTI-corrected exposure time 
taken from the EXPTIME header keyword. Finally, 
the merging of the spectra was accomplished with 
the task addspec. We ensured that each spec- 
trum contained the proper BACKSCAL-keyword and 
was linked to its corresponding background-file. 
The source spectra were individually grouped with 
dmgroup until a signal to noise ratio (S/N) of at 
least 3 per channel was achieved. A ll spectra were 
fitted with XSPEC (jArnaud I [l996h ver. 12.3.1, 



using an APE C collisional ioniza tion equilibrium 



(CIE) model ([Smith et al. 1120011 ) 



3. Results 

3.1. Imaging Analysis 

3.1.1. Counts Images 

In Fig.[TJ we show counts images of NGC604 for 
the different energy bands on a linear min/max- 
scale. These images are not exposure corrected 
and therefore preserve integer counts to allow an 
evaluation of the number of counts in each pixel. 
Pixels with observed counts consistent with back- 
ground alone are shown in white (95% Poission 
confidence). That is, assuming Poission statistics, 
we are 95% confident that pixels with < 2 counts 
(panels a and b) and < 1 counts (panels c and d) 
belong to the background. 

From these images it is obvious that NGC604 
emits significant diffuse X-ray emission, most of 
it in the soft energy band, roughly forming the 
shape of an 'U'. There is also significant soft 
emission in the north at RA = 0l' l 34 m 33.30 s , 
Dec = +30°47'42'.'0 (all coordinates in this pa- 
per are given for J20 00.0), which was previ- 
ously not detected by M aiz-Apellaniz et al. 
(J2004J). Emission seen in the medium X-ray 



band largely coincides with that seen in the soft 
band, but does not reach the same extent. In 
the hard X-ray band, we detect only two signifi- 
cant emission features, one at RA = 0l' l 34 m 35.00 s , 
Dec = +30°47'12"0 (15 net cts) and the other 
at RA = 0l' l 34 m 34.00 s , Dec = -|-30 o 47'00"0 (6 net 
cts). The bright extended source visible in the 
broad band at RA = 0l' 1 34 m 36.50 s and Dec = 
+30°47'15'/0 (ChASeM33 sour ce 346) is most 



cated at RA = 01 /l 34 m 33 s and Dec =30°46'40"0 
can be identified with SNR GKL98-94 C D'Odorico 
et al Il980l: iGordon et al. Ill998l: Ghavamian et 
al. 120051) and i s listed as source number 342 in 



likely a background AGN (see iPlucinskv et al 
2008, for details). The southernmost source lo- 



Plucinskv et alj (|2008l ). Another feature is the 
region of reduced X-ray emission running roughly 
north -south through the center of NGC604 and 
can best be seen in panels a and d. 

3.1.2. Contour Maps 

In Figs. -Ef, we provide X-ray contour maps 
for NGC604 based on exposure corrected images 
in the same energy bands as shown in Fig. [1] 
These maps were smoothed in DS9 using a ker- 
nel radius of 2 pixels to enhance the contrast be- 
tween source and background and are overlaid on 
the HST Ha image to allow a direct compari- 
son with structures of the warm ionized gas. For 
each energy band, six contour levels are plotted 
on a sqrt-scale, starting at 3tr above the back- 
ground. The b-band contours cover 4.55 x 10~ 8 
to 1.67 x 10~ 7 counts/s/cm 2 , the h-band covers 
1.68 x 10~ 8 -2.72 x 10~ 8 counts/s/cm 2 , the m- 
band covers 6.69 x 10~ 9 -3.3 x 10 -8 counts/s/cm 2 , 
and the s-band covers 2.41 x 10~ 8 to 1.67 x 
10~ 7 counts/s/cm 2 . 

The broad-band emission extends ^250pc 
along the N-S axis and ~ 240 pc along the E-W 
axis, respectively. As can be seen from the HST 
Ha image, the two hard X-ray sources as well 
as the X-ray-bright spot at medium energies do 
not have Ha-emitting counterparts. The lack of 
X-rays at the center of NGC604 is likely due to in- 
ternal absorption (see below). SNR GKL98-94 is 
not detected at energies above 2.6 keV. However, 
to identify all those features and to constrain the 
origin of the gas, it is necessary to discuss the 
Chandra observations in the framework of other 
multi-wavelength data. 

3.1.3. Multi-wavelength Comparison 

The multicolor image in Fig. [3^ combines nar- 
rowband HS T imaging (Ha and [Nil]) with the 
soft X-ray emission detected by Chandra. It un- 
ambiguously demonstrates that all bubbles and 
cavities are filled with hot X-ray emitting gas. 
The hot gas appears to be nd around the Ha ridge 
(|Tenorio-Tagle et aljliooch . roughly forming a 'U'- 
shaped structure (cf. Fig. [5b). Fig. [3] also shows 
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Fig. 1. — Counts images of NGC604 (binned by 4x4 pixels) for the broad (b), hard (h), medium (m), and 
soft (s) energy band shown on a linear gray scale. The integer counts in each band are plotted for the 
minimum and maximum count values neglecting the AGN source visible eastward of NGC604 in panels a-c. 
White represents pixel ranges belonging to the background (95% Poission confidence) . In case of NGC604 
10" corresponds to 36 pc. 
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Fig. 2. — Panel a): Solid lines encircle regions with X-ray emission from which spectra were extracted. The 
dotted line in the inlay marks the location of the Ha-'cap'. Panel b): Sk etch of NGC604, wit h regions 
and cavities defined in the text. The locations of stars and stellar clusters ( Hunter et al. Ifl996l ) and the 
'U'-shaped morphology of the X-ray gas (white dashed line) are also indicated. Panels c)-f): Contour 
maps for the different energy bands (see text). In all panels, except panel b, the underlying image is the 
continuum-subtracted HST Ha image. 
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Fig. 3. — Panel a): A multi-wavelength comparison between Ha (red), [Nil] (green), and soft X-rays 
(0.35-1.1 keV, blue) reveals that the bubbles visible in Ha are filled with hot gas and host a massive stellar 
population. It also shows the negative morphological correlation between the different gaseous phases. Panel 
b): This three color high resolution (2") composite image provides the first detailed look of the X-ray emission 
and its spatial extent. Red represents the soft, green the medium, and blue the hard band, respectively. 
Panel c): The stellar continuum emission detected in the FUV indicates the existence of young stars and 
further strengthens the link between hot gas and stars as ionizing sources. Panel d): Three of the Ha-bright 
bubbles can also be seen in the NIR at 8/im. The emission seen in this band traces molecular dust grains 
and implies that the interior of these structures has bee n excavated by stellar activity . Yellow circles are 
X-ray sources taken from the ChASeM33 source catalog (jPlucinsky et al.ll2008l ). 
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the negative morphological and spatial correlation 
between soft X-rays and Ha-emitting gas. A sim- 
ilar negative correlation is also seen in compar- 
ing the soft X-ray emission with the 8.4 GHz radio 
continuum map presented by Churchwcll & Goss 
(|1999l ). This kind of morphological anticorre- 
lation agrees well with our understanding of the 
assembly and composition of a multi-phase ISM. 
A much more detailed view of the complex mor- 
phology of NGC604 is provided by the first high- 
resolution (2") three color X-ray image of NGC604 
(Fig. [3p). On this image the 'U'-shaped soft X- 
ray emission is clearly visible, bending around the 
prominent ridge of ionized hydrogen and molecu- 
lar gas (Fig. EH). 

As pointed out by iTenorio- Tagle et alj (|20Qdh . 
the ionization structure (and thus the radiation 
field) as well as the velocity field in the western 
part is completely different from that in the east- 
ern part. This can be attributed to the presence 
of the Ha ridge which may be shielding the east 
from the violent stellar evolutionary processes and 
the hard radiation field in the west. This kind of 
dichotomy is also reflected by the distribution of 
the FUV emission (Fig. [3}:) which is a tracer of 
the stellar continuum emitted by hot stars. The 
bulk FUV emission originates westwards of the Ha 
ridge and c oincides with the nort hern part of star 



cluster A ( Hunter et al. 19961 ) which contains 



about 200 O/WR stars. Considering the stellar 
winds of such stars, it appears to be reasonable 
to attribute the existence of the X-ray emission in 
that region to stellar mass loss. However, the X- 
ray emission in the north and east of NGC604 is 
unlikely to be explained by the same kind of ion- 
ization mechanism, simply because of the absence 
of similarly massive clusters. 

As there is no significant hot gas detectable be- 
yond the outer extent of the warm ionized gas, one 
might think of NGC604 as a giant boiling Hn- 
region whose bubbles partly broke out, merged, 
and formed the observed blister-like structures. 
The released X-ray gas, however, still has not yet 
penetrated the outer H n-shell. 

3.1. 4-- Individual Regions 

Some of the X-ray emitting regions found by 
our survey o verlap with areas investig ated in 
the optical by ITenorio- Tagle et alJ ll200()h. Mafz- 
Apellaniz et al. (|2004l ). and iHunter et al 



( 19961 ). Because these studies did not use a com- 
mon nomenclature and the borders of the inves- 
tigated areas are unknown in some cases, and 
because our study also incorporates new regions, 
we decided to define our own regions and nomen- 
clature. The interesting areas are highlighted by 
solid white lines in Fig. [2^, and encircle the emis- 
sion detected in the broad-contour map (Fig. [2fc). 
Shapes and sizes, i.e. the assumed geometries, of 
the emitting regions are shown as white dashed 
lines and were estimated by examining the avail- 
able multi-wavelength data. Obvious examples are 
bubbles B1-B3 whose Ha morphology implies a 
roughly spherically-symmetric geometry. We also 
investigate some non-spherical regions which we 
call cavities (C1-C3). These cavities are also im- 
portant, because the bubbles seem to be either 
connected to or embedded in these regions, sug- 
gesting that the cavities are also influenced by the 
stellar wind and radiation field. A nice example 
is Bl which seems to be part of the much larger 
cavity CI (see Fig. [2p for a sketch of NGC604). 
Cavities C2 and C3, finally, are of interest because 
they contain X-ray emission which apparently be- 
longs to neither of the other regions. 

Regions Bl and B2 appear to be coincident with 
re gions 1 and 2 investigated by Tenorio- Tagle et 
al. d2onnh. The same areas are labeled 'A' and 'B' 
in iMaiz-Apellaniz et al. " ( 20041 ). Moreover, re- 
gion ' C and parts of region 'D' in Mafz-Apellaniz 
et al. (|2004f ) correspond to our regions B3 and 
CI , respectively. Most of the stars investigated 



by IHunter et al. 



(|l996l . their clusters A and B) 
are coincident with our regions Bl, B2, and CI. 
Cluster A covers regions Bl and B2 while cluster 
B is coincident with region CI. Regions C2 and 
C3 have no corresponding counterparts and were 
newly introduced by us. In the following we briefly 
comment on the most important observational re- 
sults for each of these regions. 

Cavity CI: From Fig. [3^ it appears that this 
region represents the largest coherent structure in 
NGC604 formed by stellar activity. CI hosts bub- 
ble Bl which contains the bulk of cluster A. In 
addition, it also harbors cluster B which contains 
about 8 09-stars (see Fig.[2jD for its location). As 
the multi-color X-ray image of Fig. [3] implies the 
whole cavity is filled with soft X-ray emitting gas. 
The noticeable depression at RA = 0l' l 34 m 32.0 s , 
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Dec = +30°47'10'.'0 is apparently caused by ab- 
sorption along the Bl shell. 

The diffuse X-ray emission is significantly en- 
hanced westward of the Ha ridge, nicely following 
that structure. This feature is also detected at 
energies above 1.1 keV and indicates significantly 
higher plasma temperatures. 

Cavity C2: Soft X-ray emission is also clearly 
detected in C2, a region in the eastern part of 
NGC604, which separates Bl and B2 from B3. No 
emission is visible on the X- ray map presented by 
Maiz-Apellaniz et al. I (|2004l ) . The western part of 
C2 contains a significant fraction of the stellar con- 
tent of cluster A. This area also emits at medium 
X-rays, as can be seen in Fig. [51 From the multi- 
wavelength data we identify 3 stellar clusters in 
the eastern part of C2 (see Fig. [2Jd for their posi- 
tions). 

Cavity C3: Another striking feature in NGC604 
visible in Fig. [2f is the soft diffuse X-ray emission 
found in the northernmost region C3 (cf. also 
Fig. [3Jd). This emission was previously not de- 
tect ed at a significant level bv Maiz-Apellaniz et 
al. 12004) and appears to be confined towards 
the north by a 'cap'-like structure seen in Ha (see 
Figs. [2k and[3b). The full extent of the X-ray and 
Ha emission can best be seen from Fig. [2J;. 

Another feature detected at medium energies 
is located north of the soft X-ray source at RA 
= 0l' , 34 m 33.5 s , Dec = +30°47'44"0. From our 
multi-wavelength data no point sources are vis- 
ible which could be responsible for the detected 
X-rays. It should be pointed out that the X-ray 
gas in C3 fills up the ga p in the radio contin uum 
distri bution observed bv lChurchwell fc Goss I (see 
19991 ). 



Bubble Bl: Bubble Bl is adjacent to the 
largest cavity CI and contains a substantial frac- 
tion of cluster A. A multi- wavelength comparison 
between FUV, X-ray, optical, and NIR emission 
(Fig. 03 shows that ChASeM33 source No. 339 
is coincident with the stars located in Bl. These 
cluster stars are also a source of strong FUV emis- 
sion and can be associated with extended thermal 
IR emission, which also hints at the presence of a 
young stellar population. In addition, this region 



also emits strongly in soft and medium X-rays. 

Bubble B2: Bubble B2 is located south of Bl at 
the southern end of NGC604, right above the SNR 
GKL98-94 (source 342 in Fig. [3]). Its morphologi- 
cal appearance is very similar to bubble Bl. This 
region also contains stars from cluster A and is 
clearly detected in soft and medium X-rays. 

Bubble B3: Bubble B3, the clearest-defined 
bubble, is visible in Ha as a luminous limb- 
brightened closed shell on the western side of 
NGC604 and is filled with soft X-ray emitting 
gas. In this case the observed morphology sug- 
gests that the warm ionized gas surrounds the 
HIM and possibly prevents it from breaking out 
into the ambient medium. No shell fragmenta- 
tion is seen, as might be caused for example by 
high internal gas pressure or thermal conduction 
due to direct exposition of the warm gas to the 
HIM. It also appears unlikely that outflow has yet 
occurred, because the motion of the gas should fol- 
low the density gradient of the ambient medium, 
that is towards the northeast and away from the 
Ha wall in the west. Since there are no signs of 
a breakout, it is justified to consider the HIM to 
be bounded by the cooler Ha gas. The other two 
bubbles, especially Bl, seem to have released their 
hot gas into cavities C1-C3. 

Optical data show bubble B3 to be marginally 
expanding, wit h expansion velocit i es of about 
Wexp < 5km/s (|Tenorio-Tagle et al.llioool ). This 
suggests that the expansion of the bubble pro- 
duces only a moderately strong forward shock as 
it sweeps up the ambient ISM. This is also sup- 
ported by enhanced [Sii]/Ha and [O ill] /Ha line 
ra tios along the outer edge o f the shell (see Fig. 1 



Tenorio-Tagle et al.ll200oh . All other line ratios 



for B3 are low (<0.33) and imply low excitation. 
Electron densitie s are likely < 100 cm" 3 f Maiz- 
Apellaniz et al. 120041 ) and therefore within the 
low density limit. Fig.[3j3 also shows that B3 emits 
significant X-rays at medium energies from a ring- 
like structure. If this structure represents the 
limb-brightened edges of a bubble, B3 also seems 
to possess (similar to Bl and B2) an additional 
temperature component which is hotter than the 
one traced by the soft emission. A higher tem- 
perature would also be consistent with our above 
assumption that the hot gas is still trapped in B3. 
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Table 1: Best fit parameters from the multi-tcmpcraturc CIE model 
ID # of dof nhp I K x K 2 i W x 



ID #of xld dof nhp 
cnts 

(1) (2) (3) (4) (5) 



CI 
C2 
C3 
Bl 
B2 
B3 
T 



645 
452 
162 
109 
126 
229 
1961 



0.723 
0.801 
0.247 
0.479 
0.560 
0.437 
1.032 



16 
18 
10 
6 
7 
15 
49 



0.77 
0.70 
0.99 
0.82 
0.79 
0.97 
0.41 



K x K 2 
lO^cm" 5 
(6) 



4 190+ - 469 

1-i»U_Q gg 4 

L573 +0.570 



0.687 H 



0.432 



-0.378 
+0.265 
0.178 
9 ooi +0.558 

1 1R7+0 085 
1.10' -0.796 

S45+ - 463 

U.O4O_ 320 
8 5Q7+3-284 
°- 0ij ' -2.634 



1-0.257 
; -0.186 

0.5801^22 
206+ ' 280 

U.ZUD_ .094 

245+° 103 

U.Z4O_ 234 

ili +a231 

u - iii -0.104 

n ior+0150 

U.1»U_q 078 

1 3Q4+ - 316 
1 -°^-0.264 



n OO+0.04 
"•■"-0.03 
n 9 7 +0.05 
u - z ' -0.04 

27 +al8 

u - z ' -0.25 



kT 2 
keV 
(7) 



kT 



+0.17 

-0.10 
+0.26 



0.18 
0.22 
0.32 
0.29 



+0.53 
-0.21 
+0.08 
-0.06 
+0.16 
-0.10 
+0.03 
-0.06 



0.80 

1-01-0.04 

0.69^1? 

1+0.11 

-0.23 
+0.23 
-0.13 

Q8+ - 23 

u - y °-0.28 

QQ+ - 20 
U - Ja -0.19 



0.61 
0.70 



u - w -0.05 
-,+0.13 
-0.04 
!+0.23 
-0.21 
2+0.39 
-0.22 
3+0.11 
-0.08 

U - 04 -O.16 
^+0.08 
-0.10 



0.55] 
0.44j 
0.28^ 
0.33j 



0.49j 



Notes: All spectra were fitted with XSPEC, using a photoelectric absorber (phabs) and a 2-temperature CIE model 
(apec, ISmith et al. I l|200ll)) and fixing the Hi column density Nn to an average value of 1.1 x 10 21 cm" 2 ( ChurchweU 
& Goss 1 19991 : iLebouteiller et al. II2006T) . In the following, subscripts refer to the lower and higher temperature 
component. Col. (1): ID of the spectral extractions regions shown in Fig. [2^,. All regions are point source-removed 
and trace the pure diffuse emission. Col. (2): The number of counts within the extracted region. Cols. (3)-(5) list 
the reduced % 2 \ the degree of freedom (dof), and the null hypothesis probability (nhp) of the best fit. Col. (6): 
Normalization constants of the fit. Col. (7): kT is the mass- weighted average plasma temperature. 



3.2. Spectral Analysis 

Despite an exposure time of about 300 ks, the 
total counts within the individual regions are rel- 
atively low (see Tabled]). We therefore grouped 
the spectra shown in Figs. 2] and [5] to attain a 
S/N > 3. All spectra are fitted with XSPEC assum- 
ing a photoelectric absorber (phabs), to correct 
for absorption along the line of sight, and a ther- 
mal plasma model. The underlying assumption 
of a thermal gas is substantiated by the fact that 
none of the alternative models (which consisted 
of a photoelectric absorber plus up to two power- 
law models) could even roughly approximate the 
observed spectrum (reduced x 2 > 3-8). 

For the thermal comp onent, we started w ith a 
single CIE (APEC model, ISmith et al.~ll200ll ). but 



this resulted in statistically unacceptable fits with 
reduced x 2 > 2.5. The quality of the fit could be 
substantially improved by adding a second APEC 
model (reduced x 2 <1.1). In addition, we ex- 
perimented al so with a non-equilibriu m ionization 
(NEI) model ijBorkowski et al.ll200ih . as the un- 
derlying physics of a fast adiabatically expanding 
and continuously shocked wind might also hint at 
NEI conditions of the gas. Compared to the two 
component CIE model, the NEI-fits produced in 
all cases consistent temperatures and very similar 
X 2 values, differing by less than 10%. However, 
since the HIM is unlikely to be represented by a 



constant temperature, as the NEI model assumes, 
but rather by a temperature distribution, we con- 
sidered the multi-temperature CIE model to best 
represent the thermal plasma. 

The Hi column density was set to TVh = 
l.lxl0 21 cur 2 (see IChurchwell fc Gossl Il999l : 
20061 ) and not allowed to 



Leboutciller et al. 



vary during the fitting process. Relative element 
abundances for NGC604 were take n for He, N, 
O, an d S from IDiaz et al. 1 (Il987fl and Vilchez 
et al. (|l988l ) 



and the Ne abundance value was 



adopted from Crockett et al. ( 20061 ). All other 



abundances were left at their revised solar values 
(|Asphmd et al. 1120061 ) and scaled by Z = 0.65Z Q , 
the value of the relative oxygen abundance. The 
best-fit statistics for each region are provided 
in Table Q] and a representative 90%-confidence 
contour plot for the (kTi,kT 2 ) parameter com- 
bination is given in Fig. [5] for the total diffuse 
emission of NGC604. The average plasma tem- 
perature kT was calculated by weighting the two 
temperature components kTi (* = 1, 2) with a fac- 
tor of \[Kil(\JK\ + \JK 2 ), i.e., an approximate 
mass-weighting of the temperatures. The average 
plasma density n c is simply the sum n\ + n 2 . 

As can be seen from the total diffuse spec- 
trum (Fig.EJ), the diffuse X-ray gas can be reason- 
ably well fitted with a thermal two-temperature 
plasma model (reduced x 2 = 1-03) which provides 
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Region B1 phabs*(apec+apec) 



Region B2 phabs*(apec+apec) 




> 2x10-' 
£ 1.5x10-3 

=5 
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" 10" ; 
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M 

]1 5x1CH 



Energy (keV) 
Region B3 phabs*(apec+apec) 





6x10- 3 



4x10- 



£ 2x10- 



1 

Energy (keV) 
Region C1 phabs*(apec+apec) 




Fig. 4- 
2.0 keV. 



Energy (keV) Energy (keV) 

Spectra extracted from regions highlighted in Fig. They cover the energy range from 0.35- 



strong evidence that the bubbles arc filled with 
coronal gas. The total, absorption-corrected, dif- 
fuse X-ray luminosity of NGC604 (0.35- 2.5 keV) 
is L x = 1-43+0^6 x lQ 36 ergs _1 . However, this is 
only 14 % of the total diffuse X-ray luminosity of 
30 Dor ( Townslev et al. Il2006h and emphasizes the 
difference between the two largest H n-regions in 
the Local Group. 

4. Discussion 

4.1. The Origin of the X-ray Emission 

There is striking evidence that NGC604 con- 



type stellar population (e.g. 


, D'Odorico & Rosa 


1981: Drissen et al. 1993: 


Hunter et al. 19961: 


Gonzalez-Delgado & Perez 


20001: Bruhweiler et 


al. 2003). It is also generally accented that 



NGC604 is a giant stellar wind-blow n mass loss 
bubble fe.gjRosa fc D'Odorico Ill982l: Y ang et al. 
ll996l : lGonzalez-Delgado fe Perez ll200oT ) powered 



by these stars. If this is true, the young stellar 
population should be responsible for the X-ray gas 
mass and luminosity. In addition, the X-ray emis- 
sion should be spatially correlated with stellar as- 
sociations as the thermal plasma is considered to 
be pr oduced by a sh ocked stellar wind (e.g. Weaver 
et al. ll977l : IShull fc Saken Ill995l ). 

With this work we provide an important test 
of the 'stellar wind'-hypothesis by investigating 
which of the regions outlined in Fig. [S^t are cur- 
rently being powered by mass loss from stars. The 
spectral fits allow us to constrain densities and fill- 
ing factors of the HIM and to calculate gas masses 
and luminosities for the individual regions. By 
constraining the mass loss history of NGC604, we 
estimate the stripped-off gas mass and check if it 
is consistent with the observed X-ray mass. Fi- 
nally, we derive the lumi nosity expected fro m the 



standard bubble model dChu et al~ _ 19951: Shull 



Saken Il995l ) and investigate whether it agrees 
with our observations. 
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Region C2 phabs*(apec+apec) 



Region C3 phabs*(apec+apec) 



: 



3x10 3 - 




1 2 
Energy (keV) 

Total: phabs*(apec+apec) 




0.5 1 2 

Energy (keV) 



_g1.5x10- 3 - 




1 2 
Energy (keV) 



Total: Confidence Contours (90%) 




<o r 

° 0.2 0.25 0.3 

Parameter: kT 1 (keV) 



Fig. 5. — Spectra extracted from regions shown in Fig. [^t cover the range from 0.35- 2.0 keV, except the 
one extracted from the total region as it also contains emission beyond 2keV. 90%-confidence contours of 
(fcTi, kT2) are plotted for the total diffuse emission. 



In cases of Bl, B2, and C2 we find good correla- 
tions between regions containing massive O-stars 
and X-ray bright regions. Given the Ha and dif- 
fuse X-ray distribution, Bl, B2, and CI could be 
connected, implying that these bubbles were pro- 
duced by the combined stellar winds from O and 
WR-stars sweeping up the ambient medium. In 
this regard, the X-ray enhancement along the Ha 
ridge in CI could be produced by a reverse shock in 
the wind after impinging on t he Ha wall. However, 
Gonzalez-Delgado fc PereTI (|2000f ) argue that the 
mechanical wind luminosity is sufficient to form 
bubbles Bl and B2 but is too low to account for 
the creation of CI. We therefore speculate that CI 
is an older structure formed by stellar winds from 
a previous star forming event; additional heating 
may have been provided by SNe. 

Furthermore, as blowout from regions Bl, B2, 
and CI is confirmed by optical spectroscopy 
(|Tenorio-Tagle et all l2000h the X -ray emission 



seen in cavity C3 could be due to hot gas stream- 
ing from Bl and B2 through CI up into C3. Such 
motions are expected to be supersonic and should 
produce some characteristic emission when the 
X-ray gas interacts with the cooler Hi and Hll 
gas. The emission feature at medium energies in 
C3 could well represent one of these signatures. 
Like the X-ray enhancement near the Ha ridge, 
this harder emission could be caused by a reverse 
shock after the upstreaming hot plasma hits the 
Ha-'cap' which seems to confine NGC604 towards 
the north. 

What kind of source powers bubble B3 and 
produces the X-ray emission: SNRs or stellar 
winds from O or WR-type stars? Perhaps the 
puzzling off-centered X-ray feature detected in 
the hard energy band at RA = 01 /l 34 m 35 s and 
Dec = +30°47'13'.'93 can provide some hints. This 
emission is significant (S/N = 3), is well aligned 
with the shell, and coincides with the low excita- 
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tion region observed in the optical (Te norio-Tagle 
et alj2000; iMaiz-Apellaniz et al. 2004). There are 
net 12cts in the hard band (see also Fig. [IJ which 
are centered at about 3.7 keV. Given that the Pois- 
son probability of having 12 or more background 
counts in the same region is -C 1%, it is highly un- 
likely that this feature is caused by a background 
fluctuation. The reality of this feature is further 
substantiated by the fact that it is visible in two 
of our four ObsIDs. Unfortunately, the counting 
statistics are too low to obtain a decent spectrum 
and to reliably constrain its nature. The multi- 
wavelength data shown in Fig. [3] and the 8.4 GHz 
radio conti nuum map p resented bv Churchwell & 
Goss (Il999h do not reveal an obvious source which 
can be associated with this emission. Although a 
contribution from SNRs cannot be ruled out, the 
fact that the emission is relatively hard and peaks 
at ^3.7keV makes these sources rather unlikely 
candidates. 

If this emission originates from NGC604 its lu- 
minosity should be in the range of Lx = 0.9 — 
1.2 x 10 34 ergs _1 (using the Portable Interactive 
Multi-Mission Simulator (PIMMS v3 . 9df] and as- 
suming power law photon indices > 1.0). This 
luminosity is consistent with the one expected 
from an X-ray binary, a pulsar wind nebula, or 
an AGN. If this source is not part of NGC604 it 
could be an AGN (seen through the shell of B3) 
with an unabsorbed flux in the 2.0-8.0 keV energy 
band of F s» 2.5x 10~ 15 ergs" 1 cm~ 2 . This flux is 
based on a peak line o f sight H I colu mn density of 
N R = 1.5xl0 21 cm- 2 (|Newtonlll980l ) and matches 
that of a typical AGN detected by the Chandra 
deep field surveys. Considering that NGC604 cov- 
ers ~ 4 x 10~ 4 square degrees on the sky, 0.84 AGN 
are expected to cover this area ( Brandt & Hasinger 
120051 ) which corroborates the AGN hypothesis. If 
we include the other AGN candidate w 17" east- 

the area in- 



2008), 



ward of B3 (jPlucinskv et al. 
creases to~8.8xl0~ 4 square degrees and 1.8 AGN 
are expected in this field. 

The 'stellar wind'-hypothesis remains to be dis- 
cussed for region B3. The multi-color image in 
Fig. shows a small off-centered stellar associ- 
ation and several other stars close to the shell 
of B3, though projection effects may need to be 
considered. The location of the association at 



the SE edge of the shell (see '*' symbol in B3 
of Fig. is also coincident with the detected 
FUV emission feature in the archival GALEX data 
shown in Fig. [3b. The GALEX FUV filter has its 
peak transmission at about 1500A, which traces 
stellar emission down to an effective temperature 
of ~2xl0 4 K. Such a temperature roughly cor- 
responds to stars with spectral type B0 or ear- 
lier. Moreover, the existe nce of young hot stars 
i n B3 is also confirmed by M aiz-Apellaniz et al 



(|2004l ) who suggest that B3 contains several late 
OB stars. However, it appears questionable that 
B3 is the result of current stellar wind activity, 
simply because the number of stars seems to be too 
low. The current energy input of these stars makes 
breakout of hot X-ray gas into the surrounding 
medium very unlikely. Therefore, we conclude 
that the energy which is required to create this 
bubble and produce the hot gas was mainly pro- 
vided by previous stellar winds or by SNe. We 
will investigate B3 and the other regions in detail 
in the following sections and test whether stellar 
mass loss is sufficient to account for the observed 
X-ray luminosity and gas mass. We will also de- 
termine whether the required number of stars is 
consistent with observations. 

4.1.1. The X-ray Gas Mass 

The determination of the X-ray gas mass al- 
lows us to investigate whether NGC604 is indeed 
powered by ~ 200 O/WR-type stars whose heavy 
mass loss drives a steady wind into the ambient 
ISM. For that purpose, we assume that the derived 
X-ray luminosities listed in Table [5] are entirely 
produced by a reverse shock hitting the wind and 
heating the gas up to X-ray temperatures. This is 
the emission mechanism assumed by the standard 
model of stellar wind-blown bubbles (Castor et al. 
ll975l : IWeaver et al. Hl977l : IShull fc Sakenlll99.^ . 

The determination of the X-ray gas mass also 
requires knowledge of the emitting volume V, the 
X-ray filling factor /x, and the electron density n e . 
These variables are related to the normalization 
constant K of the spectral fit as follows: 



K = 



10- 



4nDl 



n e riH dV, 



(1) 



4 htt p : / /exc . harvard .edu/toolkit/pimms.jsp 



where Da is the angular size distance to the source 
in (cm), and n e and njj are the electron and hydro- 
gen densities per cm -3 , respectively. The ionized 
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Table 2: Derived physical parameters based on fits presented in Table Q] 



ID 



(1) 



T 
MK 

(2) 



F- 



10- 



X,abs (unabs) 
14 „„„ „-l 



L 



erg s 
(3) 



cm 



X,abs (unabs) 

10 35 erg s _1 
(4) 



V 


M x 


n e 


(cm 


' 3 ) 


10 60 cm 3 


M Q 




fx 




(5) 


(6) 


0.1 


0.5 


0.8 


15.5±7.75 


2760±1040 


0.65 


0.29 


0.23 


5.56±2.78 


1160±490 


0.77 


0.34 


0.27 


4.62±2.31 


610±220 


0.46 


0.20 


0.17 


1.40±0.28 


630±190 


1.65 


0.74 


0.58 


0.88±0.18 


320±140 


1.34 


0.60 


0.47 


3.17±0.63 


590±260 


0.68 


0.31 


0.24 


31.1±8.58 


5780±1670 


0.67 


0.30 


0.24 



CI 
C2 
C3 
Bl 
B2 
B3 
T 



6-38lS:S 
5 10+ 2 - 67 

d.lU_ 2 44 
o 9c -+4.52 
°-^°-2.55 

O 0O + 1.28 

°- oo -0.93 
,+2.09 
-1.86 
j+0.93 
'-1.16 



6.26] 

5.6Sj 



1.96 (2.98) 
1.15 (1.61) 
0.38 (0.52) 
0.30 (0.42) 
0.46 (0.72) 
0.59 (0.82) 
4.86 (7.16) 



+0.30 
-0.40 
+0.25 
-0.27 
+0.13 
-0.14 
+0.21 
-0.26 
+0.16 
-0.28 
+0.18 
-0.26 
+0.51 
-0.28 



3.90 (5.94) 
2.29 (3.21) 
0.76 (1.04) 
0.60 (0.84) 
0.92 (1.43) 
1.18 (1.63) 
9.68 (14.3) 



+0.60 
-0.80 
+0.50 
-0.54 
+0.26 
-0.28 
+0.42 
-0.52 
+0.32 
-0.56 
+0.36 
-0.52 
+ 1.02 
-0.56 



Notes: Col. (1): ID of the spectral extraction regions shown in Fig. [2Jl Col. (2): Mass- weighted average plasma 
temperatures. We converted from keV to K using T (K) = 1.16 x 10 7 fcT (keV), with kT being the energy of a 
thermal source (see Table [TJ. Cols (3) and (4): Absorbed and unabsorbed fluxes and luminosities, listed for the 
0.35 -2.5 keV energy band. All uncertainties are given on a 95%-confidence level using the error command in XSPEC. 
X-ray luminosities assume a distance to M33 of D = 817 kpc dFreedman et al" 1 1200 lT) . Cols (5) -(6): Estimated 
volume and X-ray gas mass. Mx is calculated from electron densities if a filling factor of fx = 0.8 is adopted. 



gas mass traced by X-rays can be written as: 
M x = 1.15 m H n V fx, 



(2) 



where 1.15 x mn is the mass per hydrogen atom, 
accounting for the contribution of He and assum- 
ing a fully i onized gas w ith solar abundances (As- 
plund et al. l2006h . 

As these quantities cannot directly be mea- 
sured, reasonable assumptions need to be made. 
We are aware that these assumptions can be a se- 
rious source of uncertainty, especially in the case of 
V where most of the uncertainty is due to poorly 
constrained viewing angles, non-spherical geome- 
tries, projection effects, and (most importantly) 
unknown linear dimensions along the line of sight. 
However, despite these uncertainties the following 
calculations provide valuable insights into the en- 
ergetics and the evolution of NGC604. Moreover, 
most of the above mentioned uncertainties can be 
minimized, e.g., by a visual inspection of the avail- 
able multi-wavelength data and by providing rea- 
sonable uncertainty ranges. 

In the following we assume that bubbles Bl- 
B3, the clearest defined structures in NGC604, 
can be approximated by spheres (see Fig. [5^) with 
radiH of rei = 22.5 pc, rs2 = 19.3 pc, and rB3 = 
29.6 pc, respectively. If the bubbles are spheri- 



cally symmetric, the uncertainties associated with 
the radius would be negligible (<4%). However, 
to account for possible asymmetries and for un- 
certainties in the distance to M33, we assume an 
uncertainty in the volume of 20%. 

Estimating the volumes of the cavities CI - C3 
is much harder because of their seemingly irregu- 
lar shapes. However, if we consider that NGC604 
is e mbedded into a large reservo ir of cold gas 



Maiz-Apellaniz et al. 20041; Tosaki et al 



120071 ) through which star formation propagates 
and carves out the blister- like cavities, one can 
approximate such regions by more regular geome- 
tries. After an inspection of the HST Ha image 
we decided to approximate cavity CI by a box 
of about 140 pc x 64 pc x 64 pc in size. In order 
to correct for possible projection effects and the 
unknown 'real' extent along the line of sight, we 
adopt a 50% uncertainty in the volume (which will 
also be the default value for C2 and C3). Because 
region Bl seems to be part of CI, we subtracted 
the volume of Bl from the volume of CI. Cav- 
ity C2 can be approximated by a box which mea- 
sures ~47pc x 86 pc x 47 pc, whereas C3 is repre- 
sented by a box whose extent amounts to about 
47 pc x 58 pc x 58 pc. All resulting volumes are 



5 The dimensions were determined from the HST WFPC2 



Ha image, which provides a spatial resolution of ~0"1 
(=0.36pc). 
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listed together with their uncertainties in Table [2] 
We now proceed to constrain the filling factor 
of the hot phase and derive the electron density 
within the individual regions. Unfortunately, both 
quantities are not known a priori, but as follows 
from Eq. (2), n c ~ f x 1 ^ 2 , which allows us to calcu- 
late the density for different filling factors. A high 
value of fx is implied for the bubbles, because the 
X-ray plasma seems to be confined by the cooler 
gas (cf. Fig. [3k). Furthermore, a filling factor 
of 0.1 was derived by iGonzalez-Delgado fc Perez^ 



(|2000h for the warm ionized gas in NGC604. This 
leaves ~10% for the warm neutral and molecular 
medium; a reasonable value if one takes into ac- 
count that the stellar wind and radiation field ion- 
ized and swept away most of the neutral gas within 
these bubbles. This argument is confirmed by the 
IRAC 8/xm image shown in Fig. [3ji, in which all 
bubbles appear to be relatively free of dust. A high 
value of the X-ray filling factor fx is also expected 
for the cavities, because the shocked wind-driven 
gas should be able to easily fill a large volume (in 
the present evolutionary state cooling is negligible, 
see below). Therefore, it appears reasonable to as- 
sume a high filling factor of fx — 0.8 for NGC604 
(densities for fx — 1.0 are up to 0.04cm -3 lower). 

The resulting electron densities and gas masses 
are given in Table [21 Mass-related uncertainties 
are calculated by propagating the errors for vol- 
ume, filling factor, and density, taking into ac- 
count that Am ~ y/AKi/(AfxAV). The un- 
certainty associated with the filling factor was set 
to A/ x = ±0.1. To calculate M x ,t, the total X- 
ray emitting mass in NGC604, we simply used 
the sum of the individual volumes. The mass in- 
tegrated over the individual regions is consistent 
with the total mass within the calculated uncer- 
tainty (£M x ,id = 6100 ± 1200 M Q vs. M x ,t = 
5800 ± 1700 M©). 

4-1-2. The Stellar Mass Loss History 

In this section we investigate whether the ob- 
served stellar population can account for the X- 
ray emission by its mass loss history and which 
regions are likely to be powered by these sources. 
This assumes that the X-ray emission is produced 
by a reverse shock hitting the stellar wind. An- 
swering this question requires detailed knowledge 
about the number and type of resident stars, their 



age, and their mass loss history. 

Photometric observations of NGC604 that were 
conducted with the HST revealed that about 186 
stars are located in bubbles Bl and B2. The ma- 
jority of the se stars are of type 09.5 or younger 
cluster A in Hunter et al. Ill996l; Bruhweiler et al. 
2003). These stars are considered to have an age 
of 3 M yr and formed during an instantan eous star- 
burst ( Gonzalez-Delgado fc Perez"! [2000). Among 
these 186 star s, about 14 are classi fied as spectral 
typ e WR/Of dDrissen et al. Ill993l; Bruhweiler et 
al. l2003j). For simplicity, we assume constant 
mass loss and neglect evolutionary effects as the 
stars move off the main sequence. 

A typical mass loss rate of an 09.5 star with 
an initial mass of 40 AL i s M = 1.1 x 10 -6 Af Q /yr 



(|Lamers fc Leitherer I Il993h and of a WN7 
star with th e same initial mas s is M = 3.8 x 
10~ 5 M o /yr (jNugis fc Lamers II2OO0I ). We also 
made the assumption that the average time 
spent in the WR-phase am ounts to ~2.58xl0 5 yr 



(jMaeder fc Mevnet 1119941 ). With these values, we 



obtain ^750 Mq of material that has been released 
by the stars. This value needs to be corrected by 
the significant amount of gas which flows into the 
cavity from the shell boundaries due to evapora- 
tion via thermal conduction. In order to estimate 
this effect we make the simplifying assumptions 
of a constant plane-parallel flow, neglecting mag- 
netic suppression (kq = 1) and turbulent mixing, 
and assume that the evaporated gas is immedi- 
ately heated up to X-ray temperatures. 

Following Shull & Saken (1995), the mass gain 
rate due to thermal evaporation can be written as: 



:U. ^2^i? b CT b 5/2 



«0, 



(3) 



where /i = 0.62 tor is the mass per particle assum- 
ing a fully ionized plasma with 10% helium and 
K(T) = C T^ 2 is the thermal conductivity with 
C = 6 x 10~ 7 ergs -1 cm -1 K~ 7 / 2 . The tempera- 
ture of the bubble can be approximated as: 



5.3xl0 6 L^ 35 n 2/35 t- 6 / 3 V /7 K, (4) 



and i?b, the radius of the bubble, is given by: 

R h -65.9 Lli 5 no 1/5 1 3 6 /5 pc. (5) 

Here t 6 represents the time in units of 1 Myr and 
n is the ambient gas density per cm 3 . The 
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mechanical luminosity of the wind (in units of 
10 38 ergs" 1 ) can be expressed by: 



^38 



0.01267V* M 6 vl 



2000' 



(6) 



where iV* is the number of O-type stars having 
initial masses > 40M©, Mq is the mass loss rate in 
units of lO _6 M0/yr and U2000 represents the ter- 
minal wind velocity in units of 2000 km/s. Conse- 
quently, the mechanical luminosity produced by 
the 09 and WR-stars in Bl and B2 amounts 
to na = 2.38 (vonnn = 1.0. Lamers & Lei- 
therer Il993h an d Las .wR = 3.10 (U2000 = 0.68, 



Nugis fe Lamers 1120001 1 . respectively, which totals 



L 38 = 5.48. 

Setting n = 1.0 cm" 3 (|Rosa fc Solf 1 1 1984 ), 
Kq = 1, and te — 3 yields i?b = 179 pc and 
Tb = 6.5 x 10 6 K. These values are inconsistent 
with the observed ones and imply an evaporated 
mass of about M e — 25, 000 Mq in 3 Myr, a value 
which is more than 20 times larger than the es- 
timated X-ray gas mass in Bl and B2. Mafz- 
Apellaniz et al. (|2004l ) also found a mismatch be- 
tween observationally derived quantities and pre- 
dictions from the standard bubble model. They 
argued that the external pressure from the ambi- 
ent ISM and outflows of hot gas could slow down 
the expansion of the bubbles. 

In order to resolve the observed inconsisten- 
cies we adopt for Bl and B2 a mass-weighted 
temperature of Tb = 3.45 x 10 6 K and a bub- 
ble radius of i?b = 26.5 pc (based on Vbi+B2 = 
2.28 x 10 60 cm -3 and a spherical geometry). Un- 
der these assumptions, the evaporated gas mass 
amounts to ~770 Mq which, together with stellar 
mass loss, yields ~1500M Q in 3 Myr. This is still 
much higher than the total observed gas mass of 
~95O±25OM . 

Because we rule out significant contributions 
from SNe (see Sect. 4.2) stellar mass loss re- 
mains the favored mechanism to account for the 
X-ray gas. One way to reconcile these results is 
to consider the possibility that some of the X-ray 
emission in cavities CI and C3 is from gas orig- 
inating from the stellar wind of clusters A and 
B which streamed upwards and now fills the re- 
maining cavities. In order to investigate this, we 
now have to inc lude the additional 8 09-type stars 
from cluster B ( Hunter et al. 19961 ) in our calcu- 
lations. Among the 194 stars, a bout 15 are clas 
sified as spectral type WR/Of (jDrissen et al. 



19931 ; iBruhweiler et al. 1120031 ). The 179 O-stars 
alone lose approximately 590 Mq of their mass 
within 3 Myr, whereas the 15 WR-stars, most of 
which a ppear to be type WN7/WN8 ( D'Odorico 



& Rosa 198lh . release 150 M© into the ambi- 



ent medium. Because WR-stars evolve from O- 
stars, we need to estimate the mass loss from 
the progenitor star, too. Within 2.74 Myr this 
amounts to about 45 Mq. The 194 O/WR stars 
should have lost about 800 Mq during their life- 
time which together with an evaporated gas mass 
of M e = 3200M Q totals ~4000M Q . This mass is in 
good agreement (within 7%) with the integrated 
mass of 4300 ± 1100 M Q in regions Bl, B2, CI, 
and C3. It should be noted that the evaporated 
gas mass is a conservative upper limit, because the 
Spitzcr conductivity does not take into account 
saturation effects and magnetic fields which can 
significantly suppress conduction. 

A third important mass loading effect, pho- 
toevaporation from molecular clouds, can be ne- 
glected. Although a large molecular gas reservoir 



has been reported for NGC6 04 (e.g., iBlitz Ill985l : 



Maiz-Apellaniz et al. 1 120041 ) . radiative cooling in 



bubbles with radii r > 30 pc is expected to be 
unim portant in view of th e low density of the hot 
gas (|Shull fc Sakenlll995l ). 

We think it is possible in principle to get match- 
ing X-ray gas masses if cavity C2 is also included. 
However, there are no indications that Bl and B2 
are connected to C2. It is conceivable that the 
X-ray emission seen in C2 is in part produced by 
a stellar wind emanating from B2 caused by some 
stars of cluster A and in part by three smaller star 
clusters (see Fig.[2jD for their positions). Moreover, 
the faint FUV emission seen in C2 (in Fig. [3J:) is 
also indicative of stellar continuum and supports 
the existence of young massive stars. In the west- 
ern part of NGC604 only ~20% of the gas mass 
stems from mass loss. If we assume an even lower 
fraction of 5% for C2, 60Mq need to be produced 
by stellar mass loss, which would require 18 O-type 
stars (assuming t$ = 3 and Mq = 1.1). However, 
these stars are not observed in the HST images of 
this region, and so the gas is unlikely to be pro- 
duced by the current generation of stars there. It 
is possible that the gas is left over gas from previ- 
ous stellar winds and that SNe re-energized the gas 
and provide additional heating over time. High 
resolution spectrophotometric data are needed to 
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investigate the stellar populations and its age in 
the eastern part of NGC604 to estimate their con- 
tribution to the overall X-ray emission. 

As also noted bv lMaiz-Apellaniz et al. ( 2004 ), 
the case of bubble B3 is different from Bl or B2 
because this region might be shielded by the Ha 
ridge from the stellar wind of the central clus- 
ter. Interestingly, this bubble is expanding at 
very slow rate, indicating that this structure is 
much older than, e.g., regions w estwards of the Ha 
ridg e dTenorio-Tagle et al]|2 000: Mafz-Apcllaniz et 
al. I2OQ4J). Adopting again the same age and 
mass loss rate as above and assuming that 5% 
(~30M Q ) of the X-ray emission from B3 are pro- 
duced by a shocked stellar wind, 9 09-type stars 
would be needed to account for the observed X-ray 
gas mass. Considering that this bubble contains 
only a small unresolved stellar association and ~3 
isolated stars, presumably O-stars or red super- 
giants, the number of stars can hardly account for 
the X-ray gas mass. We therefore speculate that 
the hot gas in bubble B3 stems from a previous 
stellar wind and/or was heated by SNe. Because 
SNRs e xpanding into a low densit y gas are hard to 
detect dMcCrav fc Kafatosll 19871; Chu & MacLow 
1990), the lack of detected SNRs in this region 
is not surprising. As a result of the expected low 
energy input and marginal expansion velocity of 
B3, we also expect the bubble to remain stable 
and confine the hot plasma in its interior. 

It should be noted that there is a significant 
spread in the reported value of M for the O /WR- 
stars, depending on the adopted model. For 
our calculati ons we assumed a simpl e homoge- 
neous w ind (lLamers fc Leithererl 19931: Maeder & 
Meynet 199 41') rather than a clumped one f Nugis 



& Lamers |2000| ; iFullerton et al. 1120061 ). How 



ever, even in the case of a clumped wind in which 
the mass loss rate for O stars can be one order 
of magnitude lower, consistent gas masses are ob- 
tained within the uncertainties. 

4-1.3. X-ray Luminosities 

If the hot gas in NGC604 is produced by a 
shocked stellar wind, the observed X-ray luminosi- 
ties should also be in agreement with theoretical 
predictions from corresponding models (e.g., Shull 
& Saken Il99.4 IChu et al. Ill995h . We therefore 
calculate the expected X-ray luminosities for each 
region starting with the main volume which con- 



sists of Bl, B2, CI, and C3. We assume again 
that the observed X-ray emission is the result of 
the shocked stellar wind. 

The X-ray luminosity of a bubble can be writ- 
ten as: 

L X = (8.2 x lO^ergs- 1 )^ 077 ^ 777 ^ 1677 /^)^ 

(7) 

where v e is the expansion velocity of the bubble, 

£ is the metal abundance in solar units, and r 

is the ratio of the minimum temperature and the 

central temperature of the bubble. 7(r) is a di- 

mensionless integral which can be approximated 

5.5 



by: J( T ) 
with t = 



3.79 - 5rV 2 



0.16 £ 



-8/35 
37 1 



-2/35,6/35 

'0 r 6 



(5/3)r 3 - (5/lljr 
(see[Ch u et al 



19951 . for details). 

With L 37 = 58 it follows that r = 0.084 and 
hence 7(r) = 2.34. These values together with 
a spectroscopically determined expansion velocity 
of v 36kms~ 1 C slits Cb)-fe) in Rosa k. Solf 
Il984l : lYang et al. I Il996l ) yield a~predicted X- 
ray luminosity of Lx — 8.70 x 10 35 ergs _1 . This 
value is well within the uncertainty limit of the 
unabsorbed X-ray luminosity integrated over the 
western part of L x = 9.25±?;^ x 10 35 ergs^ 1 . If a 
clumped stellar wind is assumed the corresponding 
luminosity amounts to Lx — 8.42 x 10 35 ergs _1 . 
This finding together with a matching gas mass 
fully supports the wind-blown bubble scenario for 
the western part of NGC604 and does not require 
additional heating from SNe. 

In a similar way we derived luminosities for 
the slowly expanding regions B3 and C2 (v e 



lOkms" 1 , iTenorio-Tagle et all l2000h . They 



6.3 x 10 33 ergs" 1 and L X ,C2 
respectively. These values un- 



amount to Lx,B3 
1.1 x lO^ergs- 1 
derestimate the observed luminosities by about a 
factor of 30 and clearly rule out current stellar 
mass loss as the main contributor of the X-ray 
emission in these regions. It also corroborates our 
previous claim that the stellar populations in C2 
and B3 are not powerful enough to account for 
the determined gas mass. Although mass loss of 
the current stellar generation could account for 
approximately 4% of the observed X-ray luminos- 
ity in C2 and B3, it appears more likely that most 
of the mass of the hot gas was produced by previ- 
ous stellar wind activity with possible additional 
heating by SNe. 

It is important to note that we also checked 
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how the X-ray luminosities and gas masses ob- 
tained for the western and eastern part of NGC604 
change if the A/r is fixed to a lower (Galactic 

( Dickey fc Lockman 



1990 ) and to a n upper value of 1.3xl0 21 cm 2 



(INewton lll980h . In both cases the results were 



consistent with those obtained for the fixed av- 
erage Nh of l.lx 10 21 cm~ 2 . The maximum de- 
viations are 14%. In addition we also allowed 
the Nh to vary during the fits. The fitted val- 
ues of the Nn were poorly constrained with large 
uncertainties and a large spread in the derived 
values. For four regions (C2, B2, B3, and T) 
the fitted values of the Nh were consistent with 
l.lxlO 21 cm~ 2 . For regions CI and C3, the fit- 
ted values (0.24xl0 21 cm~ 2 and O.llx 10 21 cm" 2 ) 
were less than the Galactic column and therefore 
considered unreliable. For region Bl, the fitted 
value was 3.83x 10 21 cm~ 2 , which is larger than 
the total column through M33 and therefore also 
unreliable. Given the poor constraints on Nn from 
the X-ray spectral fits, we decided to adopt the re- 
sults of the spectral fits which had assumed a fixed 
value of l.lx 10 21 cm -2 . However, the luminosities 
and gas masses for an unconstrained iVu were al- 
ways within the range of uncertainty provided by 
the fits which assumed a lower and upper TVh of 
5.21xl0 20 cm -2 and 1.3x 10 21 cm~ 2 , respectively. 

4.2. Contributions from SNe 

We need to constrain the number of SNe to 
evaluate whether they contribute significantly to 
the energy injected into the ISM and to the de- 
termined X-ray gas mass. This kind of investiga- 
tion can only be done reliably if high resolution 
spectrophotometric data for all regions are avail- 
able. Unfortunately, the available data restrict a 
detailed exploration of the stellar population ex- 
clusively to t he western part of NGC604. Bruh- 
weiler et al. (120031) investigated about 50 OB- 
stars in cluster A covering the mass range from 
20 < M Q < 150. They assumed a Salpeter IMF 
with a — —2.3 and compared the predicted num- 
ber of stars with their observations. In the mass 
range between 80 < M Q < 120 two stars would 
have been expected, but none were detected. We 
therefore assume that these stars , if formed, ex- 
ploded as SNe during the last 3 Myr. Interestingly, 
the 3 most massive (120 < M© < 150) and ex- 
tremely rare stars still seem to exist which would 



imply that they formed during the last 3 Myr. 
However, these authors do not rule out that these 
stars are unresolved massive binaries. 

Because the above mentioned study covers only 
~25% of the massive stellar population in cluster 
A, we need to scale the expected number of SNe 
accordingly. Since the number of exploded stars 
of the remaining population is unknown and may 
also vary from region to region, we simply extrapo- 
late linearly and obtain a conservative upper limit 
of 8 SNe during the last 3 Myr. 

The total thermal energy of the gas in the west- 
ern part provided by stellar wind activity amounts 
to about 2.5 x 10 52 erg, which is roughly equivalent 
to the energy input of about 25 SNe. This num- 
ber of SNe is three times higher than our estimate 
of the expected number of SNe and is therefore 
inconsistent with the assumption that SNe con- 
tributed significantly to the heating of the gas. 
Besides the fact that the X-ray luminosity can 
be well explained by stellar mass loss alone, there 
is another aspect which argues against contribu- 
tions from S Ne. The s tarburs t models from Lei- 
therer et al. <| 19921 ) show that the total energy 
input from stars formed during an instantaneous 
starburst 3 Myr ago, is entirely provided by stel- 
lar winds. This lets us assume that no or only a 
few SNe occurred during the lifetime of the cur- 
rent stellar generation and that mass load by SNe 
is negligible. Although there could have been a 
couple of SNe, they would not have provided sig- 
nificant heating. 

The thermal energy of the X-ray gas from stel- 
lar winds amounts to ^5.0 x 10 51 erg and is five 
times smaller than the total thermal energy. This 
discrepancy could be explained, e.g., by radiative 
cooling due to excessive mas s injection by evapo- 
ration (jShull fc Saken 1995 ) or by energy losses 
from the hot interior due to outflowing gas. Be- 
cause the stellar population in the eastern part is 
completely unconstrained, we refrain from making 
predictions of the expected number of SNe. 

4.3. The overall picture 

In X-rays as well as at optical wavelengths a 
clear dichotomy is seen in NGC604. The west- 
ern part seems to be a young (~ 3 Myr) and dy- 
namically active region in which ongoing stellar 
mass loss is most likely responsible for the ob- 
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served X-ray gas mass and luminosity. This part of 
NGC604 is X-ray brig ht compared to the bub bles 
of the LMC studie d bvlChu fc MacLowl <jl990h and 
Wang k Helfand I |l99ll)- Such X-ray 3 right bub- 
bles usually require additional energy input from 
SNRs in order to reconcile the observed luminosity 
with the predicted one from the standard bubble 
model. Hence, the western hemisphere of NGC604 
seems to be a new case in which an X-ray-bright 
bubble is produced by stellar mass loss and the 
current stellar population either has not formed 
SNe, or any SNe so far were inefficient in heating 
the gas. 

The eastern hemisphere appears to be an aged 
structure in which the gas seems to be in a dynam- 
ically relaxed condition. This relaxation could be 
caused by the lack of ionizing sources and the pres- 
ence of the Ha ridge which most likely shields this 
region from the east. The likely reason why this 
part of NGC604 is also X-ray bright and has by far 
the highest temperatures is that radiative cooling 
is neg ligible i n regions of this size f Shull & Saken 
Il99a iTenorio-Tagle et al.ll2000f ). It also implies 
that the gas is indeed confined to C2 and B3. Con- 
trary to the west, the observed X-ray luminosity in 
C2 and B3 cannot be explained by current stellar 
mass loss alone. In this case, off-center SNRs could 
produce the additional emission fChu & MacLow 

Ii99nh . 

Due to the lack of identified young clusters, the 
X-ray emission in C3 most likely does not origi- 
nate from local stell ar clusters; clusters C and D 
( Hunter et al. lll~996l ) are too far off-center to con- 
tribute significantly. Therefore, the emission in 
C3 could originate from the shocked stellar wind 
expelled by clusters A and B. If this is correct it 
also follows that the ISM is highly porous, because 
the overpressured hot gas can stream freely some 
200 pc from Bl through CI into C3 and fill large 
cavities. Such supersonic motions can explain the 
large filling factors of the HIM and are a charac- 
teristic feature in models of wind-blown bubbles. 

These streaming motions should also cause 
characteristic shock features when the supersonic 
wind interacts with the gas of the shell wall. The 
medium X-ray emission along the Ha ridge and at 
the Ha-'cap' seen in the contour map (Fig. (2) and 
in Fig. [3] can be interpreted as fingerprints of such 
shocks. This shock scenario is consistent with 
the crossing timescale of the stellar wind which 



emerges from clusters A and B and traverses cav- 
ity CI, which is of the order of —5 x 10 5 yr, as- 
suming that the wind slows down from originally 
V2000 = 1 to a conservative value of i>2ooo = 0.1. 
A little less than 1 Myr is required for the wind 
to reach the top of cavity C3 (assuming the same 
wind speed). These timescales are much shorter 
than the cluster age, and therefore, if SNe and 
SNRs can be ignored, the wind could get contin- 
uously shocked over a period of several million 
years. 

5. Summary and conclusions 

We presented the first X-ray images of NGC604, 
that are deep enough to reveal the complex in- 
terplay among different constituents of the ISM. 
A strong negative morphological correlation was 
found between Ha and X-ray emitting gas. All 
bubbles and cavities are filled with hot coronal 
gas (3.3 — 6.4 x 10 6 K), suggesting high volume 
filling factors of about fx = 0.8 and densities 
n c < 0.6 cm~ 3 . 

A total, absorption-corrected, diffuse X-ray lu- 
minosity in the soft band (0.35- 2.5 keV) of Lx = 
1.43 x lO^ergs" 1 was found for NGC604, which 
i s about 14% of the corr esponding value for 30 Dor 
(jTownslev et al. 1 120061 ) . The X-ray gas mass in 



NGC604 amounts to about 6000 M Q , which is only 



15% of the corresponding mass in 30 Dor (jWang 



1999), and emphasizes the differences between the 
two largest Hn-regions of the Local Group even 
more. 

The sum of the detected X-ray gas mass in the 
western part of NGC604, containing regions Bl, 
B2, CI, and C3, amounts to -4300 ± 11OOM 
and the unabsorbed X-ray luminosity of Lx = 
9.251^14 x 10 35 ergs™ 1 places this region into the 
regime of X-ray bright bubbles. Both, the X-ray 
gas mass and the luminosity agree well with expec- 
tations from the standard bubble model and are 
consistent with mass loss from about 200 O/WR- 
stars. In view of an estimated cluster age of 3 Myr, 
the consistency between observed and expected lu- 
minosities, a nd the prediction s from the starburst 



models from Leitherer et al. ( 19921 ). the current 



stellar generation seems to be too young for SNe 
to have made a significant contribution to heating. 
The western hemisphere contradicts the general 
view in which the luminosity of X-ray bright bub- 
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bles is underestimated by the standard model and 
additional contributions from SNRs are required. 
The X-ray gas likely originates from the stellar O- 
type population in clusters A and B. We find that 
mass load from evaporation due to thermal con- 
duction contributes about 95% to the X-ray gas 
mass. 

In the eastern hemisphere the estimated mass 
of X-ray emitting gas amounts to ^1750±550M©. 
The sum of the predicted X-ray luminosity from 
regions B3 and C2 which one would expect from 
a shocked stellar wind during the lifetime of 
the donor stars, largely underestimates the ob- 
served unabsorbed luminosity of Lx = 4. 



,+0.62 . 
-0.75 



10 35 ergs _1 by about a factor of 30. This result 
supports the hypothesis that the luminosity of 
X-ray bright bubbles is inconsistent with the one 
predicted by the standard bubble model and ar- 
gues for additional contributions from SNRs. We 
rule out a shocked stellar wind produced by the 
current stellar generation as the main source of 
the X-ray emission. We cannot rule out that the 
gas was created during an era when the stellar 
population was powerful enough to drive a con- 
tinuous wind and/or SNe heated the gas. The 
number of expected SNe in C2 and B3, respec- 
tively remains unconstrained due to the lack of 
spectrophotometric data of the stellar population. 

The dichotomy between east and west was also 
fo und for the warm ionized gas. Tenorio-Tagle et 
al.HoOO) noted that the dynamics and ionization 
of the gas in the eastern part of NGC604 is com- 
pletely different from the western part. As in the 
case of the HIM, the Ha ridge appears to play an 
important role by separating both hemispheres. 
The eastern region seems to be an old part of 
NGC604 in which the warm gas appears to be in 
a dynamically relaxed condition. In contrast, the 
western region is dominated by young stars which 
create a highly dynamical ISM via strong stellar 
winds and contribute essentially 100% to the X- 
ray gas mass and luminosity. The freely-streaming 
winds in NGC604 should produce characteristic 
emission features when they interact with the shell 
walls. Therefore, we interpret the medium X-ray 
emission along the Ha ridge and at the Ha-'cap' 
as direct fingerprints of such shocks. Provided this 
is correct, these spots of enhanced emission would 
confirm the expected supersonic motions of the 
gas. 



This work has ma de use of SAOImage DS9 
(j Jove fcMandel 1120031 ) . developed by the Smithso- 
nian Astrophysical Observatory, the FUNTOOLS 
utilities, and the HEASARC FTOOLS package. 
We thank the anonymous referee for her/his valu- 
able comments. RT wishes to thank J. Raymond 
and H. Lee for helpful discussions and suggestions. 
RT acknowledges support under NASA Chandra 
award number GO6-7073A. TJG and PPP ac- 
knowledge support under NASA contract NAS8- 
03060. PG acknowledges support through Chan- 
dra grant GO6-7073B and PFW acknowledges 
support through G06-7073C. 

REFERENCES 

Arnaud, K. A. 1996, Astronomical Data Analysis 
Software and Systems V, eds. G. Jacoby and J. 
Barnes, ASP Conf. Series, 101, 17 

Asplund, M., Grevesse, N, & Sauval, A. J. 2006, 
CoAst, 147, 76 

Blitz, L 1985, ApJ, 296, 481 

Borkowski, K. J., Lyerly, W. J., & Reynolds, S. P. 
2001, ApJ, 548, 820 

Brandt, W. N., & Hasinger, G. 2005, ARA&A, 43, 
827 

Broos, P. S., Townsley, L. K., Getman, K., 
& Bauer, F. E. 2002, ACIS Extract, An 
ACIS Point Source Extraction Package (Uni- 
versity Park: Pennsylvania State Univ.) 
http://www.astro.psu.edu/xray/docs/TARA/ 
ae_users_guide . html 

Bruhweiler, F. C, Miskey, C. L., & Smith Neubig, 
M. 2003, AJ, 125, 3082 



Castor, J., McCray, R 
200, 107 

Chu, Y.-H.. 
510 



Weaver, R. 1975, ApJ, 
MacLow, M. M. 1990, ApJ, 365, 



Chu, Y.-H., Chang, H.-W, Su, Y.-L., & MacLow, 
M. M. 1995, ApJ, 450, 157 

Churchwell, E., & Goss, W. M. 1999, ApJ, 514, 
188 

Crockett, N. R., Garnett, D. R., Massey, P., & 
Jacoby, G. 2006, ApJ, 637, 741 



20 



Diaz, A. I., Terlevich, E., Pagcl, B. E. J., Vilchez, 
J. M., & Edmunds, M. G. 1987, MNRAS, 226, 
19 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 
28, 215 

D'Odorico, S., Dopita, M. A., & Benvenuti, P. 
1980, A&AS, 40, 67 

D'Odorico, S., & Rosa, M. 1981, ApJ, 248, 1015 

Drissen, L., Moffat, A. F. J., & Shara, M. M. 1993, 
AJ, 105, 1400 

Freedman, W. L., Madore, B. F., Gibson, B. K., 
Ferrarese, L., Kelson, D., et al. 2001, ApJ, 553, 
47 

Fullerton, A. W., Massa, D. L., & Prinja, R. K. 
2006, ApJ, 637, 1025 

Ghavamian, P., Blair, W. P., Long, K. S., Sasaki, 
M., Gaetz, T. J., & Plucinsky, P. P. 2005, AJ, 
130, 539 

Gonzalez-Delgado, R. M., & Perez, E. 2000, MN- 
RAS, 317, 64 

Gordon, S. M., Kirshner, R. P., Long, K. S.; Blair, 
W. P., Duric, N., & Smith, R. C. 1998, ApJS, 
117, 89G 

Higdon, S. J. U., Higdon, J. L., van dcr Hulst, J. 
M., & Stacey, G. J. 2003, ApJ, 592, 161 

Hunter, D., Baum, W. A., O'Neil, E. J. Jr., & 
Lynds, R. 1996, ApJ, 456, 174 

Joye, W. A., & Mandel, E. 2003, Astronomical 
Data Analysis Software and Systems XII, eds. 
H. E. Payne, R. I. Jcdrzejcwski, and R. N. 
Hook, ASP Conference Series, 295, 489 

Keel, W. C , Holberg, J. B. , & Trcuthardt, P. M. 
2004, AJ, 128, 211 

Lamcrs, H. J. G. L. M., & Leitherer, C. 1993, ApJ, 
412, 771 

Lebouteiller, V., Kunth, D., Lequcux, J., Aloisi, 
A., et al. 2006, A&A, 459, 161L 

Leitherer, C, Robert, C, & Drissen, L. 1992, ApJ, 
401, 596 



Long, K. S., Charles, P. A., Blair, W. P., & Gor- 
don, S. M. 1996, ApJ, 466, 750 

Maeder, A., & Meynet, G. 1994, A&A, 287, 803 

Mafz-Apellaniz, J., Perez, E., & Mas-Hesse, J. M. 
2004, AJ, 128, 1196 

McCray, R., & Kafatos, M. 1987, ApJ, 317, 190 

Misanovic, Z., Pietsch, W., Haberl, F., Ehlc, 
M., Hatzidimitriou, D., & Trinchicri, G. 2006, 
A&A, 448, 1247 

Newton, K. 1980, MNRAS, 190, 689 

Nugis, T., & Lamers, H. J. G. L. M. 2000, A&A, 
360, 227 

Oey, M. S. 1996, ApJ, 467, 666 

Plucinsky, P. P., Williams, B., Long K. S., Gaetz. 
T. J., et al. 2008, ApJS, 174, 366 

Rosa, M. 1980, A&A, 85L, 21 

Rosa, M., & D'Odorico, S. 1982, A&A, 108, 339 

Rosa, M., & Solf, J. 1984, A&A, 130, 29 

Shull, J. M., & Saken, J. M. 1995, ApJ, 444, 663 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & 
Raymond, J. C. 2001, ApJ, 556, L91 

Tcnorio-Tagle, G., Munoz-Tunon, C, Perez, E., 
Mafz-Apellaniz, J., & Medina- Tanco, G. 2000, 
ApJ, 541, 720 

Tcnorio-Tagle, G., Muhoz-Tuhon, C, Perez, E., 
Silich, S. & Telles, E. 2006, ApJ, 643, 186 

Tosaki, T., Miura, R., Sawada, T., Kuno, N., 
Nakanishi, K., Kohno, K., Okumura, S. K., 
Kawabe, R. 2007, ApJ, 664L, 27 

Townsley, L. K., Broos, P. S., Feigelson, E. D., 
Brandl, B. R., Chu, Y.-H., Garmire, G. P., & 
Pavlov, G. G. 2006, AJ, 131, 2140 

Vilchez, J. M., Pagel, B. E. J., Diaz, A. J., Ter- 
levich, E., & Edmunds, M. G. 1988, MNRAS, 
235, 633 

Wang, Q. D. 1999, ApJ, 510, L139 

Wang, Q. D., Hclfand, D. J. 1991, ApJ, 373, 497 



21 



Weaver, R., McCray, R., Castor, J., Shapiro, P., 
& Moore, R. 1977, ApJ, 218, 377 

Yang, H., Chu, Y.-H., Skillman, E. D., & Ter- 
levich, R. 1996, AJ, 112, 146 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



22 



